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ABSTRACT 

The four bright images of the gravitationaUy lensed quasar Q2237+0305 are being 
monitored from the ground (eg. OGLE collaboration, Apache Point Observatory) in 
the hope of observing a high magnification event (HME). Over the past three seasons 
(1997-1999) the OGLE collaboration has produced microlensing light-curves with un- 
precedented coverage. These demonstrate smooth, independent (therefore microlens- 
ing) variability between the images (Wozniak et al. 2000a, b; OGLE web page). We 
have retrospectively applied probability functions for high-magnification event param- 
eters to the observed light-curve features. We conclude that the 1999 image C peak 
was due to the source having passed outside of a cusp rather than to a caustic crossing. 
In addition, we find that the image C light-curve shows evidence for a caustic crossing 
between the 1997 and 1998 observing seasons involving the appearance of new criti- 
cal images. Our models predict that the next image C event is most likely to arrive 
500 days following the 1999 peak, but with a large uncertainty (~ 100 — 2000 days). 
Finally, given the image A light-curve derivative at the end of the 1999 observing 
season, our modelling suggests that a caustic crossing will occur between the 1999 and 
2000 observing seasons, suggesting a minimum for the image A light-curve ~ 1 — 1.5 
magnitudes fainter than the November 1999 level. 

Key words: gravitational lensing - microlensing - numerical methods. 



1 INTRODUCTION 

The QSO 22374-0305, sometimes known as Huchra's lens or 
the Einstein Cross (Huchra et al. 1985) , is perhaps the most 
remarkable gravitational lens yet discovered. It comprises a 
foreground barred Sb galaxy (z=0.0394) whose nucleus is 
surrounded by four images of a radio-faint QSO (z=1.695). 
Ground based spectroscopic observations have verified that 
all four images are similar QSO's at the same redshift (Adam 
et al. 1989). Broad band monitoring has shown that signif- 
icant microlensing events occur (eg. Irwin et al. 1989; Cor- 
rigan et al. 1991). Since the optical depth to microlensing 
is of order unity at each of the image positions (eg. Kent 
& Falco 1988; Schneider et al. 1988; Schmidt, Webster & 
Lewis 1998), the magnification effects on the source can be 
considered as a network of caustics moving across the source 
plane. Strong variation in a particular image results from the 
source either crossing a caustic or passing close to a cusp. 
Q2237-I-0305 provides a unique opportunity to study mi- 
crolensing events for two good reasons. Firstly, the relative 
closeness of the lensing galaxy and the nearly on-axis align- 
ment means that the time delay between the four quasar 



images is less than a day. Thus it is easy to separate mi- 
crolensing events from intrinsic variations of the QSO. Also, 
for this QSO the time-scale for microlensing events is typ- 
ically 30-50 days, reduced by a factor of 10 from typical 
time-scales of perhaps a year. 

Although Q2237+0305 has been monitored since its 
discovery, data collected by the OGLE collaboration over 
the last 3 years (Wozniak et al. (2000a,b) ; see also 



http: / / www.astro.princeton.edu / ~ogle /ogle2 /huchra. html ) 



has for the first time obtained data of sufficient coverage 
to clearly demonstrate smooth independent flux variation 
between the images. In this paper we provide interpreta- 
tions of several features in the image A and C light-curves, 
and discuss their implications for future microlensing. Sees. 
^ and ^ describe the published light-curves and the models 
used to interpret them. In Sec. |^ we discuss the HME classes 
associated with different light-curve features by comparing 
their parameters (such as height and maximum light-curve 
derivative) to model distributions. We then describe the im- 
plications for future HMEs and present qualitative scenarios 
in Sees. |^ and ^ 
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Figure 1. Image light-curves from the 

OGLE monitoring data. Figure from the OGLE web pag e (see 
http: / / www.astro.princcton.edu / ~ogle/ogle2 /huchra.html ) . 



2 EXISTING MONITORING DATA FOR 
Q2237+0305 

Fig. I shows the OGLE hght-curve for Q2237+0305, con- 
taining all points taken before the end of the 1999 observ- 
ing season (figure taken from the OGLE web page (see 



http:/ /www.astro.princeton.edu/~ogle/ogle2/huchra.htm]). 
Note that this data was published in Wozniak et al. 
(2000a,b). The light-curves have an unprecedented sampling 
rate (~ 1 week). The OGLE data adds to the approximately 
10 years of previously obtained, but less densely sampled 
photometry presented in Schneider et al. (1988), Kent & 
Falco (1989), Irwin et al. (1989), Corrigan et al. (1990) and 
Ostensen et al. (1996). Fig. ^ shows the complete compiled 
data set. Error bars are shown representing the published 
errors. Note that the data taken prior to OGLE is in R 
and r bands whereas the OGLE monitoring is in V band. 
The typical amplitudes of microlensing flux variation may 
therefore differ between data sets if the source size is colour 
dependent (Wambsganss & Paczynski 1991). 



3 THE MICROLENSING MODEL 

Our analysis involves the comparison of ensembles of 
microlensing models to Q2237-I-0305 light-curve features. 
Through this comparison we interpret observed HMEs and 
attempt to predict the features of future HMEs. This section 
provides a brief discussion of the models used. 

Throughout the paper we use standard notation for 
gravitational lensing. The Einstein radius of a microlens 
projected into the source plane is denoted by 770. The nor- 
malised shear is denoted by 7, and the convergence or op- 
tical depth by k. The basic model for microlensing at high 
optical depth comprises a disc of point masses (eg. Kayser, 
Refsdal & Stabell 1986) with a size such that a large frac- 
tion (> 99%) of macroimage flux is recovered (Katz, Bal- 



bus & Paczynski 1986; Lewis & Irwin 1995). To construct a 
microlensed light-curve we use the contouring technique of 
Lewis et al. (1993) and Witt (1993). For the microlensing 
models of Q2237-I-0305 used in the current work we assume 
the macro-parameters for the lensing galaxy calculated by 
Schmidt, Webster & Lewis (1998). We will use the stan- 
dard notation introduced by Yee (1988), to describe these 
images. Where required a cosmology having = 1 with 
Ho ~ 75 km sec~^ is assumed. 

We describe the microlensing rate in terms of the effec- 
tive transverse velocity, which is deflned as the transverse ve- 
locity that produces a microlensing rate from a static model 
equal to that of the observed light-curve (Wyithe, Webster & 
Turner 1999 (hereafter WWT99)). The effective transverse 
velocity therefore describes the microlensing rate due to the 
combination of the effects of a galactic transverse velocity 
and proper motion of microlenses. To calculate distributions 
of values for light-curve features associated with HMEs, we 
assume that the effective transverse velocity accurately de- 
scribes not only the distribution of light-curve derivatives 
during an HME (Wyithe, Webster & Turner 2000a), but 
also the distribution of orientations between the caustic and 
source trajectory, and hence the event duration. 

We have previously obtained the following normalised 
probability distributions. These were obtained under the as- 
sumption that the source size S <^ rjo- Evidence in favour of 
this assumption was presented in Wyithe, Webster & Turner 
(2000c). 

i) ps{S\{m),Veff), the probability that the continuum source 
diameter is between 5* and S + dS given a mean microlens 
mass (m) and an effective galactic transverse velocity v^ff 
(Wyithe, Webster, Turner & Mortlock 2000). 

ii) Pv{veff\{m)) the probability that the effective galactic 
transverse velocity is between Vef / and Vef / -I- dvef / given a 
mean microlens mass (m) (WWT99). 

Hi) Pm{{m)), the probability that the mean microlens mass 
is between (m) and (m) -f d{m) (Wyithe, Webster & Turner 
2000b). 

Pv{vef f\{jn)) and pmiim)) were computed using flat 



ip{Vtr 



dVtran), and logarithmic {p{Vtr 



Vtr 



-) 



assumptions for the Bayesian prior for galactic transverse 
velocity {Vtran)- PviveffKm)) was found to be insensitive 
to the prior assumed, however pm{{m}) showed some de- 
pendence. In the remainder of this paper we use pm{{m)) 
calculated using the assumption of a logarithmic prior. We 
note that the assumption of the flat prior raises the average 
light-curve derivative by a few percent. 

The functions ps(5|(m), Ve//), p„(we//Km)), pm{{m)) 
and the HME statistics presented in this paper were com- 
puted for different assumptions of smooth matter density, 
photometric error, and direction of the galactic transverse 
velocity. Since the probability functions referred to above 
were computed from a derivative analysis, the statistics that 
we compute in this paper are quantitatively similar for the 
different possible models. Therefore we present only results 
from models with no smooth matter, a transverse velocity 
direction along the image C-D axis and simulated photomet- 
ric errors assigned according to a Gaussian distribution with 
half widths of a = AM /2 in images A and B, and a = AM 
in images C and D. 

Both the microlensing rate due to a transverse veloc- 
ity (eg. Witt, Kaiser & Refsdal 1993), as well as the cor- 
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Figure 2. Compilation of light-curve data for Q2237+0305. 
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responding rate due to proper motions (WWTOOa) are not 
functions of the details of the microlens mass distribution, 
but rather are only dependent on the mean microlens mass. 
We therefore limit our attention to models in which all the 
microlenses have the same mass since the results obtained 
will be applicable to other models with different forms for 
the mass function. 

The determination of probability for the quantity 
v^f f (m) from the Q2237-I-0305 monitoring data is quite 
robust. However the probability for the source size is derived 
from a single poorly sampled HME. The small number of ob- 
servations describing the 1988 peak (Irwin et al. 1989; Cor- 
rigan et al 1991) introduces the potential for a systematic 
error in the source size equal to the ratio of the true event 
length and the inferred event length of ~ 52 days (twice the 
estimated rise time). This can be compared to the ~100 day 
separation of the two points that provide an upper-bound 
on the event duration. The resulting systematic error in the 
estimate of source size is therefore smaller than a factor of 
~ 2. In addition there may also be a component of sys- 
tematic error from the assumption that the 1988 peak was 
due to a single caustic crossing. The statistics presented in 
the following sections are therefore computed assuming prior 
probabilities for S assuming no systematic error, 



Vs {S\{m},Veff)dS, 

and systematic errors of x2 and x5 in S: 
p', {S'\{m},v,ff) dS' = ^ps {2S\{m},v,ff)dS 
and 

p's (^S'\{m),Veff) dS' = ips {5S\{m),Veff ) dS. 



(1) 



(2) 



(3) 



Our source size estimate was made from data collected in 
the R and r bands, while the OGLE light-curves showing the 
features that we wish to investigate are in the V band. This 
introduces the possibility for another source of systematic 
error if the source has significantly different sizes in the R/r 
and V bands. 

To investigate individual HMEs we must look at light- 
curve statistics for single images. Therefore unlike the calcu- 
lation of ps, Pv and pm, which used difference light-curves, 
intrinsic source variation may be important. This cannot be 
directly measured, however in Wyithe, Webster, Turner & 
Agol (2000) (hereafter WWTAOO) limits are placed on the 
intrinsic variability power-spectrum and it is shown that in- 
trinsic variability should not be an important consideration 
during HMEs. 



4 ANALYSIS OF THE OGLE LIGHT-CURVE 

In this section we apply simple statistics describing event 
heights and light-curve derivatives a postiori to specific fea- 
tures in the published light-curves. We also retrospectively 
apply the triggering function developed in WWTAOO to two 
observed light-curve peaks. We use the results to interpret 
features found in monitoring data. In particular we discuss 
whether the large scale variation in image C (1999) was due 
to the source having crossed a caustic or moved outside of 
a cusp. 

We note that choosing light-curve features a-postiori, 
and comparing them to a sample of models in order to draw 
conclusions regarding the type of event involved has the po- 
tential to introduce a statistical bias. However in the present 
case we feel that our method is justified because of the well 
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Table 1. Table of probabilities for model light-curve features corresponding to those in the 1999 image C event. Details are discussed in 
the text. 



Event type 


Size Error 


Total fraction 


Prob < 2 weeks 


Prob Mniax < 2 


Prob Mheight < 0.5 


Prob AAf™,i„ < 




C \y 1 




nil 


U.oo 


U.OO 


u.ou 


Cusp 


5x2 


0.27 


0.19 


0.84 


0.55 


0.50 




5x5 


0.26 


0.12 


0.87 


0.57 


0.50 




5 X 1 


0.14 


0.07 


0.01 


0.02 


0.08 


+ve caustic 


5x2 


0.18 


0.06 


0.01 


0.02 


0.08 




5x5 


0.32 


0.09 


0.10 


0.03 


0.10 




5 X 1 


0.58 


0.05 


0.04 


0.13 


0.98 


—lie caustic 


5x2 


0.55 


0.04 


0.12 


0.20 


0.98 




5x5 


0.41 


0.02 


0.45 


0.30 


0.98 



established prior knowledge that caustic crossing HMEs pro- 
duce larger amplitude, shorter period and more asymmetric 
(in time) variation than cusp related events. 



4.1 Statistics of the 1999 light-curve Peak in 
image C 

The image C light-curve shows a remarkable resolved peak 
described by ~ 35 points on separate days spanning ~ 7 
months (see Fig. ^). By December 1999 the image C light- 
curve had dropped to a level similar to the end of the 
1998 observing season. The peak is quite symmetric hav- 
ing reached a height ~ 0.5 magnitudes above the 1998 level. 
The maximum derivative reached both before and after the 
event peak was ~ 2 mags/year. Images B and D remained 
fairly constant during this period suggesting constant intrin- 
sic luminosity. The event properties are discussed in relation 
to model calculations of the probability for their values. The 
results are summarised in Tab. 0. 



4-1.1 Maximum light- curve derivative 

We have calculated the cumulative probability of observ- 
ing a maximum light-curve derivative M prior to the peak 
maximum: 

P^j{M < Mo) = 

dm J dS J dveff {ps{S\{m),Veff) Pm{{m)) 

X Pv{yeff\{m))Pj;, {M < Mo\S,{fn),v,ff)) , (4) 

using a sampling rate typical of the monitoring data (1 point 
per week). We have looked at the probability P for cusps 
{Pj^j^{M < Mo)), -ve caustic crossings {Pj^j_{M < Mo)) 

and +ve caustic crossings {P^,j^{M < Mo))- Fig. ^ shows 
these functions; solid, dotted and dashed lines correspond 
to P^i^iM < Mo), Pm_{M < Mo), and P^_^(M < Mo) re- 
spectively. The functions shown assume that the source size 
estimate S is correct. We find that a maximum derivative of 
~ 2 magnitudes per year is inconsistent at > 99% level with 
the event having been a +ve caustic crossing. In addition, 
we find that a —ve caustic crossing is excluded at the 95% 
level. A cusp cannot be excluded on the basis of the max- 
imum derivative (though it is higher than expected). Our 
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Figure 3. Plots of the probability for the maximum derivative 
preceding each a light-curve peak. Functions are given corre- 
sponding to derivatives that precede cusps (solid line) , and caustic 
crossings with disappearing (dotted line) and appearing (dashed 
line) critical images. Only peaks larger than 0.1 magnitudes above 
the local minima on both sides are considered. 



conclusions are barely changed (the —ve event confidence is 
then 90%) if we assume our source size has been underesti- 
mated by a factor of 2. If our underestimate is a factor of 5 
then only the +ve event can be excluded (at the 90% level). 
The results are summarised in column 5 of Tab. Q 



4-1-2 Height of the peak above the previous local minimum 

If the light-curve minimum preceding the 1999 image C 
event is assumed to have been at the level of the 1998 season 
then the height of the peak maximum above the previous lo- 
cal minimum is Mpeafe ~ 0.5 magnitudes. We calculate the 
cumulative probability for the difference between peak max- 
ima and the preceding local minima (Mheight)- Mheight is a 
function of sampling rate for sharp peaks and sparse sam- 
plings. We have therefore used a sampling rate of 1 point 
per week and integrated the probability over v^f / as well as 
5* and (m): 

PM(Mheight < Mpeak) = 
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Figure 4. Plots of the probability for the difference between a 
light-curve peak maximum and the preceding minimum. Func- 
tions are given corresponding to the peak heights for cusps (solid 
line), and caustic crossings with disappearing (dotted line) and 
appearing (dashed line) critical images. Only peaks larger than 
0.1 magnitudes above the local minima on both sides are consid- 
ered. 

X Pv {Veff\{m)) Pm {Mheight < Mj,eak\S, (m) , Vef f)) {5) 

As before we have 

computed the function for cusps {PMciMheight < Mpeak)), 
+ve caustic crossings {PM+{Mheight < Mp^ak)), and -ve 
caustic crossings {PM_iMheight < Mpeak))- Fig. ^displays 
the resulting curves (assuming no systematic uncertainty in 
the source size estimate). AMp^ak ~ 0.5 magnitudes is typ- 
ical if the event is due to a cusp, and is ruled out at the 
98% level if the event is +ve caustic crossing. If the event is 
a —ve caustic crossing the results are inconclusive. A very 
similar conclusion is reached if the source size has been un- 
derestimated by a factor of 2. The results are summarised 
in column 6 of Tab. 0. 

4-1.3 A postiori application of the triggering function to 
event 1999C 

WWTAOO describes a general function to determine how 
long one should wait (P) for a HME following a hypothetical 
observed light-curve derivative (T). If a light-curve deriva- 
tive Tobs^ATobs is observed, predictions F{P\Toba±A.Tobs) 
can be made about forthcoming events specifically for the 
current data using the sampling rate identical to that of the 
observations. On the 19th of June 1999 monitoring from the 
OGLE collaboration (OGLE web page) showed image C ris- 
ing at a rate of 1.21-1.78 mags/year. Since at this level P 
is not sensitive to the sampling rate for sampling spacings 
smaller than two weeks, we used observations on the 10th 
June, 19th June and 1st July. 

For the present calculation the triggering function dis- 
cussed in WWTAOO is slightly modified. Rather than search- 
ing for peaks following a derivative larger than some hy- 
pothetical value, we search for peaks following light-curve 
derivatives in the measured range. Thus since the algorithm 



Figure 5. Plots of the probability for the delay of a light-curve 
peak following the observed image C light-curve derivative on the 
19th June. Functions are given corresponding to the delays for 
cusps (solid line) , and caustic crossings with disappearing (dotted 
line) and appearing (dashed line) critical images. 

steps along the curve in steps equal to the sampling rate, cer- 
tain events (in particular -\-ve caustic crossings) are missed 
if the derivative jumps from below the lower bound to above 
the upper bound during one sampling interval. Therefore the 
results describe the relative likely-hood of observing the dif- 
ferent types of event following a current light-curve deriva- 
tive. As a consequence of the triggering algorithm searching 
for derivatives along the light-curve in one direction, then 
the model light-curve derivative at the point of the trigger 
is systematically biased toward T — AT. P is therefore an 
overestimate. This bias is minimised by reduced photometric 
error. 

Fig. Is] shows the resulting triggering functions F+ , F- 
and Fc computed assuming no systematic uncertainty in 
the source size estimate. The solid, dotted and dashed lines 
correspond to Fc{P\Tobs ± Ar<,6,), F-{P\Tobs ± ATobs), 
F+{P\Tobs ± ATobs) respectively. We find that the observed 
trigger precedes a caustic crossing HME ~ 75% of the time. 
The event peak is most likely to occur ~ 3 months following 
such a trigger if the event was a cusp, 2 months later if it 
was a ~ve caustic and 1 month later if it was a +ve caus- 
tic crossing. Note that having observed a derivative in the 
quoted range means that a +ve caustic crossing is very un- 
likely due to the rapid change in the light-curve derivative. 
The light-curve peaked at ~ 1 July, about 2-weeks after the 
observed derivative. This is surprisingly early for all types 
of events (though P is an overestimate) . At the 90% — 95% 
level, caustic crossing events are excluded. This result holds 
if a systematic source size error is assumed. The probabili- 
ties of both the class of event following the trigger as well as 
the arrival time are summarised in columns 3 and 4 of Tab. 

4.2 What class of event have we seen in image C ? 

Constraints on the event type of the OGLE image C HME 
are placed by both the maximum derivative observed prior 
to the event peak, and the height of the peak above the pre- 
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Figure 6. Plots of the probability for the difference between the 
local light-curve minima either side of the peak maximum. Func- 
tions are given corresponding to the probabilities for cusps (solid 
line), and caustic crossings with disappearing (dotted line) and 
appearing (dashed line) critical images. Only peaks larger than 
0.1 magnitudes above the local minima on both sides are consid- 
ered. 



vious minima. -|-«e caustic crossings are excluded by both 
these measurements (even when potential systematic errors 
in source size are assumed). In addition, assuming no sys- 
tematic uncertainty in source-size, —ve caustic crossings are 
excluded by the measured maximum derivative. The cusp 
interpretation is consistent with both measurements. 

Because triggers in the 1999 image C light-curves trig- 
ger were relatively small, the fraction of events predicted for 
different classes of HME by the triggering function does not 
restrict the type of event observed. In particular, the trig- 
ger was not large enough to rule out either cusps or —ve 
caustic crossings. However the OGLE data shows a previ- 
ous rise in the light-curve of image C occurring between 
~ 300 and 500 days prior to July 1999. This rise occurred 
between observing seasons, but the net change is J; 0.8 mag- 
nitudes suggesting that a - \-ve caustic crossing event may 
have been missed (see Sec. 4.4). The typical separation of 
double peaked events calculated by Witt, Kayser and Refs- 
dal (1993) is 300-500 days. This fact coupled with the com- 
parative rarity of a -f t;e event following the observed trigger 
leads to the conclusion that if the 1999 event were a caustic 
crossing then it is more likely to have been the second peak 
in a double horned event (— ue caustic crossing) than a -\-ve 
caustic crossing. However, based on the peak arrival time, a 
cusp event is more likely than a —ve caustic crossing (or a 
-\-ve caustic crossing). 

Based on this evidence we conclude that the event ob- 
served for image C was probably a cusp event rather than 
a caustic crossing. However, further monitoring may show 
that this interpretation is false, and therefore that the mi- 
crolensing models are in error. Thus the event will provide 
a valuable test of whether our current microlensing models 
can be used to reliably interpret the initial stages of HME's. 




-200 
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Time (days) 
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Figure 7. The light-curves for the 1988 image A peak (thick 
line), and the 1999 image C peak (thin-line). The two peaks have 
been placed on the same time-axis such that the peaks approxi- 
mately coincide at t = 0. 



J^.2.1 Asymmetry oj the 1999 image C event 

At the time of writing the image C light-curve was still in 
decline, but appeared to be decelerating {+ve second deriva- 
tive) . Observations of the final stages of this event will pro- 
vide a further property that can be used to distinguish be- 
tween the different possibilities. After the light-curve has 
flattened out, the difference between the light-curve min- 
ima immediately preceding and following the event can be 
measured: 

AMmin = Mmin{left) - Mmin{right). (6) 
We have calculated the distribution of these values as before: 
Pam„„. ( AAf™i„ < AMobs) = 



dm dS dv^ff {ps (S\{m),Veff) 



xp,n {{m)) Pv {Veff\{m)) 

X PAM„i„ {AMmin < AMobs\S, {m),V^ff)) . (7) 

The functions were calculated corresponding to cusps, —ve 
caustic crossings and -\-ve caustic crossings (Pc, P- and P+ 
respectively). The resulting curves are plotted in Fig. ^ (for 
the case of no systematic error in the source size estimate). 
The results are summarised in column 7 of Tab. |l[ 

While the trailing minimum has not yet been observed, 
it may provide a discriminate between the interpretations 
of the 1999 image C event as a cusp and a —ve caustic. In 
particular, if the light-curve flattens out at a level approxi- 
mately equal to that of the 1998 season {AMmin ~ 0), then 
the —ve caustic crossing interpretation will be ruled out at 
the 95% level. Assuming that the previous light-curve min- 
imum occurred during 1998, the -\-ve caustic crossing inter- 
pretation is already ruled out at > 95%. 



4.3 Comparison with the 1998 image A event 

The first unambiguous microlensing signal was the rise and 
fall of image A in 1988 during a period when the other im- 
ages remained at a consistent level (ie. no intrinsic vari- 
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Figure 9. Scatter plots of peak-height above half maximum vs. full width at half maximum. Left: peaks due to caustic crossings. Right: 
peaks due to cusp events. Light and dark dots correspond to images A and C respectively. 7^ > 0, -fc < 0. 
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Figure 8. Plots of the probability for the maximum second 
derivative on the leading side of the peak maximum. Functions are 
given corresponding to the probabilities for cusps (solid line), and 
caustic crossings with disappearing (dotted line) and appearing 
(dashed line) critical images. Only peaks larger than 0.1 magni- 
tudes above the local minima on both sides are considered. 

ability). Some insight is gained through comparison of the 
shape of this light-curve peak with the 1999 image C event. 
The two events are described by very different data sets. 
While the OGLE data has provided excellent coverage of 
the 1999 peak, allowing quantities such as the maximum 
derivative and peak height to be calculated, the 5 obser- 
vations of the 1988 peak provide a much cruder record. In 
particular the values of maximum pre-peak derivative, the 
peak height above the minimum and the difference between 
the minimum cannot be computed. On the other hand, our 
source size estimation was determined from the assumption 
that the image A peak is a caustic crossing. Statistical de- 
terminations of the event type are therefore only meaningful 
as a check of the self-consistency of the calculations. In that 
vein we note that the one leading derivative that can be 
measured from the data is T ~ 5 mags/year. Calculations 



in WWTAOO suggest that the observations should have been 
followed by a caustic crossing in ~ 2—3 weeks or a cusp event 
in ~ 5 weeks, which lie between and after the two brightest 
observations respectively. 

Fig. ^ shows light-curves of the two peaks, placed on 
the same time-axis such that the maxima approximately co- 
incide at t = 0. The image A and C peaks are shown by 
thick and thin lines respectively. The large initial rise of the 
1988 event measures the minimum to the maximum gradi- 
ent (which is surely significantly larger given the high sec- 
ond derivative). This sharp rise is not replicated in the 1999 
event which had a maximum derivative of ~ 2 magnitudes 
per year. Since we have a lower bound for the maximum 
derivative of the 1988 image A event which is consistent with 
all three types of HME, the size of the observed derivative 
cannot be used as a discriminate (although for a cusp event 
there is a 90% chance of Mmax < 5 magnitudes per year). 
However, the leading three points also measure a lower limit 
for the maximum second derivative. Fig. |^ shows probability 
Pj^ for the maximum value of second derivative on the lead- 
ing sides of cusp events (solid line) and +ve (dashed line) 
and —ve (dotted line) caustic crossings computed for image 
A (no systematic error in source size was assumed): 

Pm{M < Mobs) = 

X p„ {veff\{m)) P^j (M < Mobs\S, {m),vesf)) . (8) 

The derivatives were calculated using a sampling rate cor- 
responding to the initial 3 observations of the 1988 peak. 
Comparing the image A second-derivative of Moba ~ 0.5 
magnitudes per year per day to Fig. ^ we infer that the 
1988 peak was probably a -\-ve caustic crossing. 

A second feature to be noted from Fig. ^ is that the 
1999 peak appears to have a much longer duration than 
the 1988 peak. Fig. ^ shows scatter plots of event height 
above fuU-width-at-half-maximum (fwhm) vs. fwhm for im- 
ages A (light dots) and C (dark dots). The plot on the left 
shows the relationship for caustic crossings while the right- 
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hand plot is for cusp events. The plots highlight the rise- 
time - peak-height correlation for caustic crossings, and the 
cloud of points due to smooth light-curve variations (both 
of which were pointed out by Witt & Mao (1994)). The 
separation of the points into two categories demonstrates 
the intuitive notion that cusp events have longer durations 
than caustic crossings. The cusp event parameters shows a 
peak width lower-limit corresponding to the caustic crossing 
correlation. Therefore, the systematic bias introduced into 
the source size determination by the assumption that the 
1988 event was due to a caustic crossing can only result in 
an over estimate of source size. In addition. Fig. ^ demon- 
strates that a cusp event contains very little information on 
source size/transverse velocity due to the lack of any corre- 
lation with peak height. The inference that the 1988 peak 
was a caustic crossing and the 1999 peak was a cusp event 
is consistent with Fig. ^ 



4.4 A 1997-1998 event in the OGLE image C 
light-curve ? 

Fig. ^shows a lS.Mobs ~ 0.8 magnitude rise between the 1997 
image C minimum and the 1998 level, suggesting an event 
in between those observing seasons. We assume that the 
intrinsic source luminosity was approximately constant over 
this period, which is supported both by the facts that image 
A changed by ^ 0.2 magnitudes and that the other images 
show opposite trends. Comparing the image C change to 
the probabilities in Fig. ^ we find that AAfobs is consistent 
with a -\-ve caustic crossing having occurred between the 
1997 and 1998 observing seasons, but rules out —ve caustic 
crossings (~ 99%) and cusp events (~ 95%). We therefore 
infer that a -\-ve caustic crossing was missed between the 
1997 and 1998 observing seasons. 



5 PREDICTIONS OF FUTURE HMES 

In this section we calculate probability functions for the like- 
lihood of observing future HMEs in images A and C given 
current light-curves. 



5.1 The next image C HME 

In this section we assume that there was a ~\-ve caustic cross- 
ing between the 1997 and 1998 observing seasons, and inves- 
tigate when we should next see a caustic crossing in image 
C. These calculations follow Witt, Kayser & Refsdal (1993) 
who calculated the separations in dimensionless units of the 
different combinations of +ve and —ve events. However we 
have included both our estimates of (m) and v^ff, and the 
cusp as a third class of HME. 

Due to the typical diamond formation of fold caustics, 
the case of a +ve followed by a ~ve caustic crossing is com- 
mon. Similarly, inspection of model light-curves shows that 
cusp events follow —ve caustic crossings as the source moves 
past the cusp associated with that caustic (this feature is 
seen in the double horned profile that is characteristic of 
the Chang-Refsdal lens). However we have inferred that the 
OGLE image C light-curve shows a -\-ve caustic crossing 
followed by a cusp event. Such a combination is much less 
common and is due to the source moving past a cusp formed 
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Figure 10. Examples of double- horned profile (top) and double 
horned profile with an additional cusp event (centre and bottom). 



from two fold caustics other than the one responsible for the 
+ve caustic crossing HME. It can also be seen in model light- 
curves. Examples of the two scenarios are shown in Fig. hol 
The upper panel shows an example of a double horned event. 
The lower two panels show examples of a double horned 
event surrounding a cusp event. The source is shown pass- 
ing very close to the cusp (centre panel), partly coming in 
contact, and passing further away, producing a lower am- 
plitude event (lower panel) . The light-curves were produced 
using our most likely model for the microlensing parameters 
[veff = 300fcmsec"\ (m) = 0.1 Mq and 5" = 5 x lO" cm). 
The intervals (eg. 1-2) quoted below refer to the intervals 
between the events labelled on these plots. 

From these two scenarios, we generate the probability 
functions for 4 different event separation statistics (plotted 
in Fig. [HI ) . The left hand and right hand panels show func- 
tions computed for transverse velocities aligned with the C- 
D and A-B axes respectively. Fig. |l^ also contains a key 
corresponding to the intervals between event types shown 
in Fig. These intervals are defined below. 
Firstly, we assume that the 1999 peak was a —ve caustic 
crossing. 
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Figure 11. Solid lines: Plots of the probability for the time in months until the next caustic crossing assuming that the 1999 light-curve 
peak is a cusp related event following a +ve caustic crossing. Dashed lines: Plots of the probability for the separation in months of a +ve 
HME followed by a second HME. Dot-Dashed lines: Plots of the probability for the separation in months of a +ve HME and a second 
HME that bracket a cusp event. Dotted lines: probability for the separation of a cusp event following a —ve caustic crossing. The left 
and right hand plots correspond to transverse velocity directions aligned with the C-D and A-B axes respectively. Only peaks larger than 
0.1 magnitudes above the local minima on both sides are considered. 



oq 



Figure 12. The probability for the ratio of separations in months 
between a cusp peak (following a —ve caustic crossing and the 
nearest caustic on either side. The dark and light-lines correspond 
to transverse velocity directions aligned with the C-D and A-B 
axes respectively. Only peaks larger than 0.1 magnitudes above 
the local minima on both sides are considered. 

I) interval 1-2: The dashed hnes show cumulative proba- 
bilities for the separation of two adjacent caustic crossing 
HMEs where the first is a +ve caustic crossing. 
Secondly, the 1999 peak is interpreted as a cusp event fol- 
lowing a +ve caustic crossing. 

a) interval 1-3: The dot-dashed line shows the distribution 
of caustic separations in a double horned event where a cusp 
lies between the two caustics. In this case the typical sep- 
aration of caustics is larger due to the fact that the cusp 
is generally formed from independent caustics, so that the 
chance of a cusp lying inside a diamond or another cusp is 
higher if the separation (and therefore area) is larger. 



Hi) interval 1-4: The dotted lines represent the probability 
for the separation between a cusp and the (immediately) 
preceding +ve caustic crossing. 

iv) interval 4-3: The solid lines represent the probability 
for the separation between a cusp event (that has followed a 
+ve caustic crossing) and the subsequent caustic crossing. 
As expected Hi) and iv) are very similar. The separation of 
an inferred +ve caustic crossing between the 1997 and 1998 
seasons and the 1999 image C peak is consistent with Fig. ^ 
whether the 1999 peak is interpreted as a —ve caustic cross- 
ing or a cusp event. However if the 1999 peak is interpreted 
as a cusp, we expect another caustic crossing to follow. Fig. 
PH shows that in this case, the separation between the cusp 
peak and the subsequent caustic is likely to be ^ 5 years 
(although may be significantly longer). 

These statistics for the separation of caustic crossings 
that surround a cusp does not make use of the information 
that the observed cusp event followed the caustic crossing by 
~ 500 days. An alternative way of analysing the separation 
between the cusp event peak and the next caustic crossing 
is to calculate the ratio of the time between the cusp event 
peak and last caustic crossing and the time between the 
cusp event peak and next caustic crossing. The probability 
of this ratio is plotted in Fig. ha. The typical ratio is 1 and 
we expect it to be between ~ j and ~ 4, yielding a most 
likely arrival time for the next event of ~ 500 days, and 
an upper limit of ~ 2000 days (~ 90%). Also plotted in 
Fig. hd is the probability of the ratio for a random position 
(dashed lines). The agreement of the curves illustrates the 
independence of the positions of the cusps and unassociated 
caustics. 

5.2 The next image A HME 

In addition to the remarkable peak observed during 1999 in 
the image C light-curve, the OGLE data also shows image 
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Figure 13. Plots of the probability for the delay of a light-curve 
peak following the observed image A light-curve derivative on the 
30th October. Functions are given corresponding to the delays for 
cusps (solid line) , and caustic crossings with disappearing (dotted 
line) and appearing (dashed line) critical images. 
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Figure 14. Best fit of the near caustic approximation (Eqn. 
to the rise in the image A light-curve. 



A brightening over the entire season, with the most rapid 
variations occurring in the latter observations. We have ap- 
plied the triggering function (as described in Sec. 4.1.3) to 
the image A light-curve. 

On the 30th of October 1999 monitoring by OGLE 
(OGLE web page) showed rises in image A, of 1.41-1.88 
mags/year (we used observations on the 20th of October, 
30th of October and 9th of November). Fig. |l| shows the 
resulting triggering functions -F+, F- and Fc- Functions 
are shown assuming that the source size estimate S is 
correct. The solid, dotted and dashed lines correspond to 
Fc{P\T,bs±/^Tob,), F^{P\Tots±ATobs), F+{P\Tobs±ATobs) 
respectively. We find that the results are similar to those ob- 
tained for the June image C trigger. The observed trigger 
precedes a caustic crossing HME ~ 80% of the time. The 
event peak is most likely to occur ~ 1 — 3 months following 
such a trigger if the event was a cusp, ~ 2 months later if 
it was a —we caustic and ~ 1 month later if it was a +ve 



caustic crossing. Having observed a derivative in the quoted 
range means that a +ve caustic crossing is very unlikely, 
and a —ve caustic crossing is the most likely option. Un- 
fortunately, these results predict that an event will occur in 
image A between the 1999 and 2000 observing seasons. If 
the impending event is assumed to be a —ve caustic cross- 
ing, with a previous minimum occurring during the 1998 
season, then from Fig. ^ we predict that the image A light- 
curve should have a subsequent minimum at a level ~ 1 — 1.5 
magnitudes fainter than the November 1999 level. 

If the source is small with respect to rjo and therefore 
the inter-caustic spacing, and the brightening of image A 
is due to the imminent disappearance of a pair of critical 
images, then the rise can be modelled using the near caustic 
approximation of Chang and Refsdal (1979): The flux fp of 
a point source at a small time At = tcaust — t from a fold 
caustic is 



fp = fo + e{At)- 



(9) 



where fo is the magnification of the non-critical images, 
9{At) the Heaviside step function and ao a constant de- 
scribing the strength of the caustic. The choice of which 
points to use in a fit of this type is somewhat arbitrary. 
We have chosen the data following JD 1400, which is af- 
ter the apparent infiection in the light-curve. The fit is 
shown in Fig. ^ giving parameter values of fo — OAOmJy, 
ao = 3.50 mjy days'^ and tcaust ~ 1554 days (~ 11th Jan- 
uary 2000). The final figure is of particular value since it 
predicts the time of the caustic crossing, tcaust agrees with 
the most popular value for the caustic arrival time according 
the triggering calculation. 



6 LIGHT-CURVE INTERPRETATIONS 

While it is pointless to try and fit the observed light-curves 
directly either by a functional form (except very near a caus- 
tic crossing), or with model light-curves, it is illustrative in 
summary to draw by hand possible interpretations of the 
data. Fig. |l^ shows data for images A and C from Wozniak 
et al. (2000a,b) and the OGLE web page. The top panel 
shows possible light-curves corresponding to our most prob- 
able interpretation. The solid curve shows the inferred image 
C event at ~ 4200 days, the 1999 light-curve peak as a cusp, 
and the potential second caustic crossing discussed in Sec. 



5.1 



The dashed line shows the potential — ve caustic crossing 
expected following the rise in the light-curve at the end of 
the 1999 observing season. The lower panel in Fig. |l^ shows 
schematics of less favoured options: the interpretation of the 
1999 image C peak as a —ve caustic crossing, and the rise 
in image A heralding a +ve caustic crossing. 



7 CONCLUSION 

We have applied simple event statistics a postiori to features 
in the OGLE light-curves of Q2237-I-0305. In the specific 
case of the 1999 peak we conclude that the event was due to 
the source passing outside of, but close to a cusp rather than 
to a caustic crossing. This hypothesis may be confirmed or 
refuted when the trailing peak minimum is observed dur- 
ing the 2000 observing season. In addition, we find that the 
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Figure 15. Cartoons of our interpretations of past light-curves 
and predictions of future behaviour. 



image (j light-curve rise between tJie 1997 and 1998 UCjLt; 
observing seasons was a caustic crossing that resulted in two 
new critical images. The hypothesis of a cusp event follow- 
ing the first half of a double horned event is a rare feature. 
However its consequence is that we expect another caustic 
crossing high magnification event to follow in the image C 
light-curve. Our models predict that this caustic crossing is 
most likely to arrive in ~ 500 days, and can be expected 
within ~ 2000 days (~ 90% confidence). However it may be 
considerably longer before the next image C caustic cross- 
ing is observed, particularly if the transverse velocity moves 
perpendicular to the shear vector in image C. By applying 
the triggering function developed in WWTAOO to the rise 
in the image A OGLE light curve, we predict that a caus- 
tic crossing high magnification event will occur between the 
1999 and 2000 observing seasons which will result in the loss 
of a pair of critical images. We therefore expect that the im- 
age image A light-curve should have a minimum at a level 
~ 1 — 1.5 magnitudes fainter than November 1999. 
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